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ABSTRACT 



We determine the microlensing event rate distribution, T(tE), that is properly 
normalized by the first year MACHO group observations (Alcock et al. 1997). 
By comparing the determined T(£e) with various MF models of lens populations 
to the observed distribution, we find that stars and WDs explain just ~ 49% of 
the total event rate and ~ 37% of the observationally determined optical depth 
even including very faint stars just above hydrogen-burning limit. Additionally, 
the expected time scale distribution of events caused by these known populations 
of lenses deviates significantly from the observed distribution especially in short 
time scale region. However, if the rest of the dynamically measured mass of 
the bulge (~ 2.1 x 10 10 M & ) is composed of brown dwarfs, both expected and 
observed event rate distribution matches very well. 



Subject headings: The Galaxy — gravitational lensing — dark matter — Stars: 
low-mass, brown dwarfs 
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1. 



Introduction 



The gravitational microlensing experiments toward Large Magellanic Cloud (Ansari et 
al. 1996; Alcock et al. 1996) was originally proposed to search for Massive Astronomical 
Compact Halo Objects (MACHOs) (Paczyhski 1986; Griest 1991). On the other hand, the 
experiments toward the Galactic bulge (Alcock et al. 1997; Udalski et al. 1994; Alard 1996) 
were initiated mainly to check the feasibility of the experiments since a certain event rate is 
expected from the "known" population of lenses, i.e., stars. However, experiments toward 
the galactic bulge not only provide a new probe of the galactic structure but also become 
an important tool that can resolve the controversial "disk dark matter" problem (Bahcall 
1986; Kuijken & Gilmore 1989). 

There already have been several trials to determine the mass function (MF) of galactic 
lenses from the result of microlensing experiments (Zhao, Spergel, & Rich 1995; Zhao, 
Rich, & Spergel 1996; Jetzer 1994; Han & Gould 1996). However, these analyses either 
determined event rate distribution by assuming some arbitrary functional forms of MF, e.g., 
power law, or could provide at most very crude information about lenses. For poorly-known 
lens populations such as brown dwarfs (BDs), this type of analysis might be the only 
approach currently possible, and thus the derived MF based on the observed Einstein time 
scale distribution is very uncertain and model dependent (Mao & Paczyhski 1996; Gould 
1996). The Einstein time scale is related to the physical parameters of a lens by 



where te is the Einstein ring radius, v is the source-lens transverse speed, m is the mass of 
the lens, and D Q \, D\ s , and D os are the separations between the observer, lens, and source. 
The better approach will be, however, first estimate the expected event rate distribution 
from the well constrained lens populations, e.g., stars and white dwarfs (WDs), and then 
test other possible lens populations. By doing this, one can set the upper limits on dark 
lens populations. However, still no quantative estimate of detailed event rate and time scale 
distribution by even known population of lenses has been made. 

In this paper, we determine the microlensing event rate distribution, T(£e), that is 
properly normalized by the first year MACHO group observations (Alcock et al. 1997). 
By comparing the determined T^e) with various MF models of lens populations to the 
observed distribution, we find that stars and WDs explain just ~ 49% of the total event rate 
and ~ 37% of the observationally determined optical depth even including very faint stars 
just above hydrogen-burning limit. Additionally, the expected time scale distribution of 
events caused by these known populations of lenses deviates significantly from the observed 
distribution especially in short time scale region. However, if the rest of the dynamically 
measured mass of the bulge (~ 2.1 x 10 10 M , Zhao, Spergel, & Rich) is composed of brown 
dwarfs (BDs), both expected and observed event rate distribution matches very well. 





(1.1) 
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2. Model I: Lenses of Known Population 



2.1. 



Mass Function 



To determine the microlensing event rate distribution caused by stars, it is required 
to construct properly normalized MF of galactic stars. Unfortunately, the MF cannot be 
obtained from observation, because it is difficult to measure stellar masses. Instead, one can 
construct MF by inferring stellar masses from well constrained luminosity function (LF) by 
using mass-luminosity relation. 

The LF of galactic bulge stars is constructed as follows. For stars brighter than 
My ~ 4, most of these stars are clump giant stars, the LF is constructed based on the 
de-reddened /-band LF determined by J. Frogel (private communication). For the part of 
the LF fainter than My ~ 4, we adopt that determined by Light, Baum, & Holtzmann 
(1997) by using the Hubble Space Telescope. They find that the LF of stars in an angular 
size of arcmin 2 is well represented by 



in the magnitude range 4 < My < 10. For stars even fainter than My = 10, we extend the 
LF by adopting that of stars in the solar neighborhood (Gould, Bahcall, & Flynn 1996) 
under the assumption that galactic disk and bulge stars have similar LFs. The absolute 
V band magnitude is determined based on a distance modulus of //bulge = 14.5, which is 
equivalent to the distance of 8 kpc to the bulge. For the luminosity to mass conversion we 
use the mass-luminosity relation determined by Henry & McCarthy (1993). The relation 
has the form 



Since the mass-luminosity relation uses U-band magnitude, we convert the /-band to 
U-band magnitude by 



In addition to visible stars, we also include WDs into a member of known lens 
populations (Adams & Laughlin 1996). In our model MF, WDs are uniformly distributed 
in the mass range 0.5 M < m < 0.6 M Q and their total number is normalized so that 
there are 10 times more WDs than clump giants (stars brighter than My ^ 4) considering 
the life expectancy of clump giants and age of the Universe. 

The model LF, ^^dMj, and the corresponding MF, $ m dm, are shown in Figure 
1 and in the upper panel of Figure 2, respectively. Both functions are normalized for 
objects in a physical area of pc 2 at the galactic center, corresponding to an angular area of 
(0.43 arcmin) 2 . In the figure, the range of LF, —4 < Mj < 15, is divided into 1000 intervals, 
while that of the MF, M < m < 1.5 M , is divided into 100 intervals. 



log$ y = 0.083M y + 2.47, 



(2.1.1) 



' 0.4365 - 0.0971M y + 0.002456M 2 , 
log(m/M ) = I 1.4217- 0.1681M y , 

1.4124 - 0.2351M y + 0.005257M£, 



(M v < 10.25), 

(10.25 <M V < 12.89), 

(M v > 12.89). 



(2.1.2) 




(2.1.3) 
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Then what fraction the galactic bulge mass is composed of these stellar and WD 
lenses? These values are determined from the relation 

^- = (2.1.4) 

-^bulgc ^pop 

where Lbuigc = /bulge u dV is the total luminosity of the bulge and Ex = Jbw v d£ is the 
surface light density seen through BW. Here v is the 3-dimensional light density distribution 
of the galactic bulge and the notations /bulge u dV and J* BW vdi represent volume integral 
over the whole bulge and line integral along the line of sight toward BW, respectively. For 
the computation of m pop , we adopt Kent (1992) bulge light density distribution model of 
the form 

i \ _ J 1-04 x 10 6 (s/0.482)- L85 L G pc" 3 , (s < 938 pc), (l) 1 _v 

U[S) 1 3.53Ko (s/667) L & pc" 3 , (s > 938 pc), 1 j 

where s 4 = R 4 + (z/0.61) 4 , the coordinates (R,z) represent the galactocentric distance 
along and normal to the galactic plane, and K is a modified Bessel function. From this 
light density model, we find the surface light density of E L = 2412 L & pc~ 2 and the total 
luminosity of Lbuigc — 1-8 x 10 10 L Q . In addition, the surface number and mass densities of 
individual lens populations toward Baade's Window (BW) are computed from the MF by 

/ $ m dm, 
Jbw 

S m = / m$ m dm, (2.1.6) 

JBW 

resulting in = 3422 pc -2 and S m = 1562 M pc~ 2 for stars, and Sjv = 614 pc~ 2 and 
S m = 369 pc -2 for WDs. Then the total masses of individual populations in the bulge 
are m sta r = 1.18 x 10 10 M Q and m WD = 0.25 x 10 10 M for stars and WDs, respectively. 
According to this MF model, the combined total mass of stars and WDs in the bulge 
comprises ~ 70% of the dynamically determined bulge mass of 2.1 x 10 10 M Q (Zhao, 
Spergel, & Rich 1995). In Figure 1, we list the values of Sjv, S m , and Ex,, and resultant 

^^pop' 



2.2. Event Rate 



The next question is, then, how many events are caused by these known populations of 
lenses? The event rate distribution of bulge self-lensing events (Paczyhsky et al. 1994) for 
a single source is computed by 

Tbuigc^E) = e(f E ) J dD os n(D os ) J dD ol (2r E )n(L> ol ) J dv y dv z vf(v y ,v z ) 5(t E - t' E ) 

x JdD os n{D os ) \ (2.2.1) 
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where efe) is the detection efficiency, n(D os ) and n(D Q \) are the number densities of source 
stars and lenses, (v y ,v z ) are the components of the transverse velocity, v, parallel and 
normal to the galactic plane, and f(v y ,v z ) represents their distributions. The factor 2r E , 
i.e., Einstein ring diameter, and v are included because lenses with larger cross-sections and 
higher transverse speeds result in higher event rate. For the case of disk-bulge events, on 
the other hand, the rate distribution computation can be simplfied by 

J dD oX (2r E )n(D ol ) J dv y dv z vf(v y ,v z ) 5(t E - t' E ). (2.2.2) 

because bulge stars are located at a large distance, i.e., D os ~ 8 kpc, compared to a typical 
observer-lens separation. 

For the galactic bulge and disk matter distributions, we adopt Kent bulge model (see 
eq. [2.1.6]) and double exponential disk model. The disk matter distribution model has the 
form 

'#- 8000 z 



n(R, z) = no exp 



h R + h 



(2.2.3) 



where the values of the radial and vertical scale heights are Hr = 3.5 kpc and h z = 325 pc, 
respectively, and hq is the normalization (see below). The velocity distributions for both 
disk and bulge components are modeled by a gaussian, i.e., 



f(vi) oc exp 



2(7 



iey,z, (2.2 A) 



where the means and dispersions are (v yj( iisk, ^V.disk) = (220, 30) km s _1 and (-u^disk, o" 2 ,disk) = 
(0, 20) km s -1 for the disk component, and (v y ^ u \ gc , o^buige) = (0, 93) km s _1 and 
(v z , bulge, °z,buigc) = (0, 79) km s" 1 for the bulge component, respectively. See Han & Gould 
(1995, 1996) for more details. 

For the computation of event rate distribution it is required to know the 3-dimensional 
number density n, but what is observationally known for the bulge matter distribution 
is the light density v. Under the assumption that light is distributed the same way as 
matter is distributed, one can convert light density into number density by multiplying the 
number-to-light ratio, i.e., 

n(D ol ) = [-) v{D ol ). (2.2.5) 

We find the number-to-light ratio of the galactic bulge to be N/L = Sat/S^ = 1.66. After 
normalizing the bulge self-lensing event rate distribution by using N/ L ratio, the disk-bulge 
event rate is scaled so that the optical depth ratio between the two components becomes 

rdisk _ I[T' disk (t E )t E /e(t E )}dt E _ 0.5 x IP' 6 

Tbuigc j[T'{t E )t E /e{t E )}dt E 1.2xl0- 6 ' { ''> 



based on the optical depth computation by Han & Gould (1995). 
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Once the event rate distribution for a single star is obtained by equations (2.2.1) 
and (2.2.2), the total event rate distribution for all monitored iV* stars during the total 
monitored time T is computed by 

r(t E ) = N*TT>(t E y, V = r' bulgc + r^. (2.2.7) 

We find the total exposure of MACHO experiment by 

V2r J i=1 e{t E ,i) 

MACHO group reported 45 events toward the galactic bulge (Alcock et al. 1997). They 
determined the optical depth based on 41 events, excluding 4 events either failed the final 
cut or suspected as variable stars, to be r = 2.4 x 10~ 6 . By using the detection efficiency 
e(t E ) also provided by MACHO group, we find the total exposure to be 

N,T = 2.08 x 10 9 stars days. (2.2.9) 



The finally determined event rate distribution expected from the known lens population 
(thick solid line) is shown in the lower panel of Figure 2, and it is compared with 
the observed distribution (histogram). In the figure, the bulge-bulge and disk-bulge 
event rate distributions are represented by short and long dashed lines, respectively. 
One finds that the event rate expected from the known populations of stars and WDs, 
r exp , alone cannot explain the observed event rate distribution, r o b s - One also finds 
that the biggest difference between the two rate distributions arises in short time 
scale region. We find that events caused by stellar and WD lenses make up only 
r c x P ,tot/r ob s,tot = jr cxp (t E )dt E / JT ohs (t E )dt E = 49% of the total observed event rate, and 
they comprises r cxp /r obs = J [t E T cxp (t E ) / e(t E )} dt E / J[t E T ohs (t E ) /e(t E )a}dt E = 37% of the 
optical depth. 



3. Alternative Mass Function Models 

3.1. Model II: Additional Brown Dwarf Population Lenses 

In previous section, we show that additional population(s) of lenses are required 
to explain the observed galactic bulge event rate distribution. These candidate lens 
populations might be black holes, neutron stars (Venkatesan, Olinto, & Truran 1997), BDs, 
etc. Among these candidates, the most probable population would be BDs because the 
major difference between r obs and r cxp arises in short time scale region. 

Therefore, we make an alternative MF model in which lenses are composed of already 
known populations of stars and WDs plus BDs which make up the rest of the dynamically 
determined galactic bulge mass. For this case, BDs comprise 31% of the total galactic bulge 
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mass, corresponding to 0.65 x 10 10 M . This fraction of BDs in the Galaxy might sound 
too big according to the theory of star formation (Adams & Fatuzzo 1996; Graff & Freese 
1996). However, one cannot rule out the possibility that there are a large number of BDs 
from a different population of stars, e.g., Population III. In the MF model we assume BDs 
are uniformly distributed in the mass range 0.06 M & < m < 0.08 M Q . 

The event rate with this alternative MF model is determined similarly, but with 
different value of number-to-light ratio. This is because the surface number density increases 
compared to that of model I, while the surface light density does not change (see Table 1) 
due to the dark nature of BDs. In Figure 3, we present the model MF and corresponding 
event rate distribution. One finds that the determined event rate distribution including BD 
population matches very well with the observed distribution. The ratio between expected 
and observed event rates and optical depths are r cxPi tot/r o bs,tot = 0.96 and r cxp /r bs = 0.54, 
respectively. There are three very long time scale events with t& ^ 70 days whose nature 
is hard to understand under reasonable matter and velocity distribution models. If these 
three very long events are not included, the ratios become r exPitot /r obs tot = 1.03 and 

Texp/Vobs = 1.28. 



3.2. Model III: Power-law Mass Function 

Zhao, Spergel, & Rich (1995) claimed that the observed time scale distribution could 
be explained if the dynamically measured mass of the bulge were distributed in a Salpeter 
power-law MF between 0.08 M Q and 0.6 M , i.e., 

$ m (m) = C n m' 2 - 5 , (3.2.1) 

where C n is the normalization constant. From the surface mass density computed by 
S m = (£L/L bu ige)Mp p = 2814 M Q pc~ 2 the value of C n is obtained by 

where &' m is an arbitrarily normalized MF. The model III MF is shown in the the upper 
panel of Figure 4 (solid line) and it is compared with model II MF (dotted line). Based on 
this MF we reproduce the event rate distribution and it is shown in Figure 4 (lower panel). 
One finds that the event rate distribution matches impressively well with the observed one. 

However, this simple picture does not hold up under closer examination. First of 
all, all bulge mass cannot be in objects m < 0.6 M , since bulge MF has been measured 
for objects m > 0.6 M Q . According to our MF of stars and WDs, these surely-existing 
relatively massive objects account for important fraction of total number and mass of the 
bulge; ~ 31% of the total number and ~ 51% of the total mass of stars. Therefore, they 
should be included in optical depth and event rate computation. To make up these massive 
stars, model III MF overestimates low-mass stars, whose MF also does not match with 
observation. 
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TABLE 1 

The Lens Population Inventory 



mnHpl 


n^ra TTlPt PT* 


llTl if 


stars 


populations 
WDs BDs 


total 


T 

_L 


777 

" L pop 


1U i v -i 


1 17 

1.1 1 


0.28 




1 4 s ) 




y 


oVllPPts TiC~^ 


3422 


614 




4036 






M & pc~ 2 


1562 


ooy 




1931 




E L 


L pc~ 2 


2412 






2412 




AT 1 T 

is j Li 










1.00 


II 


'"'pop 


10 10 M Q 


1.17 


0.28 


0.65 


2.10 




v 


nViippt^ nr — ^ 


3422 


614 


10887 


14923 




/ i — 
in 


Mr,, DC ^ 


1562 


369 


871 


2802 




E L 


L pc~ 2 


2412 






2412 




N/L 










5.33 


III 


^pop 


10 10 M 


2.10 






2.10 






objects pc~ 2 


16283 






16283 




E m 


M pc- 2 


2814 






2814 




E L 


L pc" 2 


2412 






2412 




N/L 










6.75 



NOTE. — The lens population inventory of various mass function models. In the table 
m pop represents the total mass of each lens population within the galactic bulge, and Eat, 
E m , and E L are the surface number, mass, and light density seen through BW. All these 
surface densities are normalized for objects in an area of pc 2 at the galactic center. The 
number-to-light ratio is determined by N/L = Eat/E^. 
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Figure 1: The absolute /-band LF of galactic bulge stars. The LF is normalized for 
stars in a physical area of pc 2 at the galactic center, corresponding to an angular area of 
(0.43 arcmin) 2 , and its range, —4 < Mj < 15 is divided into 1000 intervals. 
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Figure 2: The MF of known population of lenses, stars+WDs (upper panel), and resulting 
event rate distribution (lower panel). The MF is normalized for objects in an area of pc 2 
at the galactic center, and its range, M Q < m < 1.5 M , is divided into 100 intervals. 
The total expected event rate distribution (thick solid line) , which is sum of bulge-bulge 
(short-dashed line) and disk-bulge (long dashed line), is compared with that of MACHO 
experiment (histogram). It is clear that the event rate distribution from the known 
population of stars and WDs alone cannot explain the observed distribution. They make 
up, respectively, just ~ 49% and ~ 37% of observed event rate and optical depth. Note the 
difference between the two distributions arises especially in short time scale region. 
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Figure 3: Model II MF (upper panel) in which lenses are composed of known population of 
stars and WDs plus BDs, which make up the rest of dynamically determined galactic bulge 
mass. The resulting event rate distribution is shown in the lower panel. The notations are 
same as Figure 2. One finds that the expected distribution matches very well with the 
observed distribution. 
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t E (days) 



Figure 4: Model III MF (upper panel) and corresponding event rate distribution 
(lower panel). The model follows a Salpeter power-law MF in the mass range 
0.08 M Q < m < 0.6 M Q . The notations are same as Figure 2. Compared to model II MF, 
model III MF does not include surely existing brighter, and thus massive, stars which make 
up nearly half of the total mass of stars. Despite this unrealistic model, it produces event 
rate distribution that matches with the observed one (lower panel). 



